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0^ ' An extension of the fermionic particle-particle propagator is pre- 

sented, that possesses similar algebraic properties to the single-particle 
Green's function. In particular, this extended two-particle Green's 
' function satisfies Dyson's equation and its self energy has the same 

analytic structure as the the self energy of the single-particle Green's 
function. For the case of a system interacting by one-particle poten- 
^f*j , tials only, the two-particle self energy takes on a particularly simple 

form, just like the common self energy does. The new two-particle self 
energy also serves as a well behaved optical potential for the elastic 
scattering of a two-particle projectile by a many-body target. Due to 
its analytic structure, the two-particle self energy avoids divergences 
that appear with effective potentials derived by other means. 
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I. INTRODUCTION 

Ph. 

Many-body Green's functions, also called propagators, provide useful tools for 
' investigating various properties of many-body systems like solids, atomic nuclei, 

C"| , atoms, and molecules. The so called single-particle Green's function is the most 

O ■ fundamental in a whole hierarchy of related propagators. One of the reasons for the 

outstanding importance of the single-particle Green's function is that it satisfies a 
well-known equation — the Dyson equation . Among the wealth of applications 
of this equation we want to mention that the Dyson equation has been exploited 
for calculations of ionisation spectra of molecules as well as scattering cross 

sections ||. In fact, the Dyson equation allows one to derive exact effective one- 
particle equations for the elastic scattering of a particle by a many-body target j],|7j . 
In such a one-particle equation, a nonlocal and in general energy-dependent po- 
tential — called optical potential — includes the correlations of the many-particle 
system. The Green's function optical potential — also called self energy — has 
several advantages upon the optical potentials derived by other means, e. g. by Fes- 
hbach projection |8|,|9| from the multiple particle Schrodinger equation. One point 
is the good-natured energy dependence that avoids divergences for high energies. 
Also the self energy closely resembles the properties of the phenomenological optical 
model potentials that have been common in nuclear physics for a long time p0| . 
Further we want to mention that the self energy of a system interacting by one- 
particle potentials takes on a very simple appearance and provides an exact optical 
potential already in first order perturbation theory. This is usually not the case 
for optical potentials derived by other means. In general, the well developed ap- 
proximation schemes for the self energy based on perturbation theory account 
for a balanced treatment of the correlation, both in the target and the scattering 
state. This balanced treatment is much more difficult to achieve in wavefunction 
approaches ||. 
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The next member in the hierarchy of Green's functions is the two-particle Green's 
function [0. This is a function of three energy variables. Although this function sat- 
isfies an integral equation with kernels that can be approximated in a diagrammatic 
approach — the Bethe-Salpeter equation — for most problems it is too complicated 
to handle. Instead, one uses propagators with only one energy variable left, e. g. the 
polarization propagator and the particle-particle propagator fis] ] which have 

been used to calculate excitation energies |l4|,[l5| and double-ionisation spectra 
|T(j| [l^ ], respectively. Furthermore, approaches have been reported starting from 
the latter propagator to calculate amplitudes for two particle scattering via optical 
potentials |ll||2(|. Although there are some analogies between these propagators 
and the single-particle Green's function, the Dyson equation is not passed on to 
these two propagators. Neither by diagrammatical analysis nor from the algebraic 
point of view is there any close analogon to this fundamental equation. 

The main goal of this paper is to show that it is possible to construct extensions of 
these two-particle propagators (only the particle-particle propagator for fcrmions is 
considered explicitly) that have the major algebraic properties in common with the 
single-particle Green's function, namely those that lead to a Dyson equation and a 
two-particle self energy. It is particularly interesting to note that this two-particle 
self energy serves as an optical potential for the elastic scattering of a two-fermion 
projectile by a correlated target. Possible applications include, e. g., the scattering 
of a Deuteron by an atomic nucleus or the scattering of a Positronium by a molecule. 
In general, the usefulness and quality of possible approximations are determined by 
the algebraic properties of the optical potential. The two outstanding features of 
our new, two-particle self energy are the good-natured energy dependence that has 
a finite limit for high energies and the fact that influences resulting from a one- 
particle potential are accounted for exactly already in the energy independent first 
order. These properties are straightforward consequences of the algebraic analogies 
to the single-particle Green's function, but they are in no way a matter of course 
with two-body effective potentials. 

This paper is organized as follows: After briefly reviewing the algebraic prop- 
erties of the single-particle Green's function we will show in section |ll] by explicit 
construction, that an extension of the particle-particle propagator is possible that 
inherits the major algebraic properties from the single-particle Green's function and 
satisfies a Dyson equation. Discussing the properties of the extended propagator 
and its self energy we also consider the special case of a system interacting only by 
one-particle potentials. In section HI we then show how the extended propagator 
can be used for deriving well behaved effective two-particle scattering equations. 
We will, however, develop these ideas in a general framework allowing for cither 
bound or free two-particle projectile states. Also, we will state explicitly the prop- 
erties of the extended two-particle propagator that are prerequisites for derivation 
of effective scattering equations. 



II. CONSTRUCTION OF TWO-PARTICLE PROPAGATORS 



In this section we will construct an extended two-particle propagator that inherits 
a couple of fundamental properties from the single-particle Green's function. Al- 
though the single-particle Green's function and the Dyson equation are well known 
and described thoroughly in many standard textbooks ]jj|T5|] , we briefly review the 
algebraic derivation of the Dyson equation |2~l| , in order to motivate our construc- 
tion of extended two-particle propagators. We carry on by applying the formalism 
developed for the single-particle Green's function to the particle-particle propa- 
gator and construct one possible example for an extension of it. This extension 
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has similar algebraic properties to the single-particle Green's function. We show 
that it is possible to derive Dyson's equation analogously and define a two-particle 
self energy. We examine the properties in relation to the single-particle Green's 
function. In particular we consider the case of a system interacting with a one 
particle potential only. In order to be able to apply the extended propagator to the 



scattering problem of section III, we Fourier transform the propagator and Dyson's 
equation into time space. 



A. Single— Particle Green's Function and Dyson Equation 

The single-particle Green's function Q in energy space is defined by: 

G pg (u>) = (*„| a p 1 ■ 4 

ui — H + E™ + irj 

<*o| aj 1 , a p |* c ) (1) 

E£ + H - vr) 



Here, \^> ) denotes the bound ground state of the correlated iV-fermion system with 
Energy . The operator a q (or a q ) creates (destroys) one fermion in an orbital 
that belongs to the one-particle quantum number q. The set of possible indices 
{q} covers the one fermion Hilbert space. The positive infinitesimal 77 provides 
the correct time ordering when the expressions are Fourier transformed into time 
space. The Hamiltonian H is the usual Fock space Hamiltonian that contains the 
kinetic energy as well as possible interactions with external potentials and two body 
interactions. As usual it is convenient to present it as a sum of a zeroth order part 
H Q and an interaction part H\. 

H = H + H X (2) 

We choose the zeroth order Hamiltonian to be a diagonal one particle operator with 
respect to the underlying orbital basis: 

H = ^eia\ai (3) 

i 

The two parts of the single-particle Green's function describe the propagation of 
quasi-particles and quasi-holes in the correlated A-fermion system. Formally we 
can interpret each of the two parts as a projection of an operator resolvent onto the 
states a p \*S? ) and aj |\& ) respectively. These states belong to the Hilbert spaces 
of N ± 1 fcrmions. They are neither orthogonal among each other nor normalized, 
since the scalar products yield one particle density matrices: 

(*„| a p a\ |* ) or (* | a\a p |* } (4) 

If, however, we formally introduce the composite states 

these form an orthonormal set with respect to a suitably defined scalar product: 1 

(Y p G \Y q G ) = a p a\ |* ) + a\a p |* > 

= (y \{a p ,a\}\* ) = 5 m (6) 
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Mathematically speaking the composite states live in the direct sum of the N + 1 
particle Hilbert space and the dual of the N — 1 particle space. 

The states \Y G ) allow the following shorthand notation for the single-particle 
Green's function: 

Gpq (uj) = {Y p G \^-\Y q G ) (7) 

The operator Ti G in the composite space is defined by: 

- G ( H - - i V \ 

H -{ EV-H + ir,) (8) 

Eq. (Q) now allows us to interpret the full single-particle Green's function as the 
projection of an operator resolvent onto a set of orthonormal states. 

These states { \ Y' )} span a subspace of the extended Hilbert space they live 
in. We want to draw attention to the point, that these states are labelled by 
the complete set of one particle quantum numbers, regardless, whether or not the 
indices refer to particle or hole states. Due to this fact, the spanned subspace has 
the same dimension as (i. e. is isomorphous to) the one particle Hilbert space. Now 
we can extend the set { |1^ G )} to an orthonormal basis { \Q G )} of the composite 
space and use matrix partitioning 2 in order to derive an expression for the inverse 
of the matrix G(lo) with elements given by Eq. (|7|). This directly yields Dyson's 
equation. We present it in the following form: 

QiwY 1 = u>l ~g~^(uj) (9) 

Here the matrix e is diagonal and contains the zeroth order single-particle energies 
e p as well as the infinitesimals ±ir]. The self energy is defined by 

E(w)=a -e + H, \j-U. ■ (10) 

■=bb 

Here H is a matrix representation of the operator Ti in the basis { \Qj)} 3 { |^ G )} 
of the extended space: 

M IJ =(Q?\n G \Qj) (n) 

By virtue of the subdivision of the basis set into two parts, the matrix is structured 
into blocks: 

Mirt:) (12) 

The block index a refers to the set { \Y G )}, the index b to the complemental part 
of the basis. The zeroth order of tL aa is given by e. For more details on the theory 
we refer to pi]] . 

It depends on the point of view, of course, whether the Dyson equation (jp) is seen 
as an equation that relates quantities which are defined separately or as a definition 
of the self energy. In the latter case, it is Eq. ([[o|) that has to be derived. We 
want to remind the reader, that historically the self energy was defined through its 
perturbation series and consequently the Dyson equation had to be proven [2^j2^| . 
In contrast to this historic background we understand the Dyson equation rather as 
an algebraic property of the single-particle Green's function. When we say that we 
try to construct two-particle propagators that satisfy a Dyson equation, we mean 
it exactly in this sense. 
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B. Two— Particle Propagators 



After this excursion to the single-particle Green's function we will focus our 
attention on two-particle propagators. We consider the particle-particle propagator 
and its constituents: 



n rs>rV H = n+ iPV (a;) - n- >pV M (13) 

KsM") = a ^ WT^FZ a °' a r' ^ 14 ) 

oj — H + hj™ + ir\ 

nr«,rvM = <*o| 44 — t^n ~ a — ( 15 ) 

The part 11+ r / s /(w) contains the relevant information for elastic two particle scat- 
tering off the correlated ground state, whereas II~ r < s /(w) can be used to calculate 
double ionization spectra. We would like to mention that when transformed to 
time-space, the particle-particle propagator is just an additive part of the full two- 
particle Green's function consisting of two specific time-orders of the latter. 

The analogy to the single-particle Green's function fails in the point that neither 
the JV ± 2 particle states 



44 |* ) or a r a s |* ) (16) 



nor the composite vectors 



44 1* ) 



(17) 



satisfy any kind of orthogonality relation. Even more, the overlap matrices do not 
have to be regular at all. This is most easily seen by considering the special case 
of a noninteracting system (zeroth order). In this case, the states defined above 
arc slater determinants. All these states, however, vanish, if the index r refers 
to an orbital that is occupied in the ground state determinant and s refers to an 
unoccupied orbital. I. e. ,the space spanned by these states is not isomorphous to 
the two particle Hilbert space but has a lower dimension. 

Comparing with the one particle case, we see, that likewise the states \^ a ) do 
not span an isomorphous space to the one particle space, but the composite states 
\Y^) do. That was the reason for being able to express the single-particle Green's 
function as a matrix inverse or resolvent in a space spanned by the complete set of 
one-particle quantum numbers. The fact that the states \Y q ) form an orthonormal 
set has the consequence that the self energy does not contain any linear contributions 
in the frequency variable u>. 

In order to restore the analogy to the single-particle Green's function we introduce 
the following extended composite states: 

/ 44 l*o) \ 

\Yrs) = Xrs (18) 

\(*o| 44/ 

With the choice 

Xrs ■= 4 l vaC > ® (*o| 4 ~ 4 l vac ) ® (*o| 4) ( 19 ) 

where |vac) denotes the Fock-space state with no particles, and the "metric" 
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(20) 

these states fulfil the following "orthogonality" relation: 3 

(Yrs | A l^r's') &rr'$ss' &rs'&sr' (21) 

This is as much as a unit matrix in a space spanned by two particle quantum 
numbers. The second term on the right hand side of Eq. ( pl[ ) is due to antisymmetry 
in the one particle indices that comes from the fermionic nature of the particles. 

The extended states live in a composite vector space, denoted by Y, which is 
the direct sum of the spaces, the three components live in. These are: the N + 2 
particle space, the dual of the N — 2 particle space, and the space Xrs lives in, which 
is the tensor product of the one particle Hilbert space and the dual of the N — 1 
particle space. We have found other possible choices for the extension Xrs that live 
in different spaces and utilize different metrics [p4[ . In this paper we only give one 
example. 

In the composite space Y, the metric (i defines a binary product by: 

(\A),\B)) = (A\fl\B) (22) 

We want to mention that due to the indefinite choice of the metric jl, this binary 
product is not a scalar product. 

In order to define a Green's function analogous to Eq. (Q), we first have to define 
an operator that takes on the role of liF from Eq. (||). A suitable choice for our 
goals is: 

H-E* -irj 

H :=\ F« +E? -H( N -V + in o j (23) 

Eg — H + i 

The symbol H^- 1 ' denotes an operator that acts like the usual Fock space Hamilto- 
nian on the first component (the one particle part) of the product space, i. e. 
H <g> I. Accordingly, H^ N ~^ acts only on the second component, which is the dual 
of the N - 1 particle space, i. e. H^ N ~^ = 1®H. 

Now we are in the position that we can define a new extended two-particle Green's 
function by 

Grs,r's'(u) ■= {Yrs |A ~ ~ \Y r >s') (24) 

uj — TL 

It follows directly from this definition that this function contains the particle- 
particle propagator as an additive component: 

H 

XV s> ) (25) 




u-HW ~E* + H( N ~ 1 ) -ir) 

= n rSjr v (w) 




4 |vac) * (* | al — r -a r |* ) 

ui — Eg + H — in 



s 



'} (26) 
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The astcrix * denotes convolution with respect to u) and the symbol {r < — ► s} 
stands for the term in braces reappearing with interchanged indices. 

The new terms that appear in the extended propagator additionally to the 
particle-particle propagator are analytic functions of w in the lower half of the 
complex plane. Right above the real axis they feature single poles and branch cuts. 
The poles appear at energies that are sums of single-particle energies (eigenvalues of 
H^) and quasi-hole energies (single ionization energies or eigenvalues of H^^ 1 ^). 
Branch cuts originate in the continuous parts of the spectra of the single-particle 
Hamiltonian and the (N — l)-particle Hamiltonian H^^ 1 ^, respectively. In 
computations with finite basis sets, of course, only discrete single poles appear. 
In general, the new poles introduced in the extension Q{lo) of the particle-particle 
propagator Tl(w) are distinct from the poles of the latter. Therefore, the quantities 
of interest, namely the poles and residues of the particle-particle propagator, can 
be identified and extracted from the extended propagator. 



C. Dyson Equation for the Extended Propagator 

The starting point for our considerations is definition (|24| ) and the orthogonality 
relation ( plf ) for the composite states \Y rs ). Since these states are antisymmetric 
with respect to a permutation of the indices r and s, we introduce the index restric- 
tion r > s in order to make the set of vectors {|y rs )} linear independent. Now we 
can extend this set to a basis { \Qi)} D {\Y rs )} of the space Y. We may choose this 
set orthogonal with respect to the binary product: 



(Qi\ A \Qj) = Sumi 

Due to the indefiniteness of the metric /t, we cannot normalize all states to ra; 
but we can require that 



mi G {-!,+!} 



(27) 
= +1, 

(28) 



We define the following matrices: 



:= bijmi = (Qi\ fi \Qj 
■■= (Qi\ m \Qj) 



(29) 
(30) 



Just like in section II A. Eq. (|12|), these matrices are structured into blocks by the 



subdivision of the basis { \Qi)} into the subset and its complement. Of 

course, the upper left block of /i is just the unit matrix due to equation (pl|): 



—a a 



i 



(31) 



After some simple manipulations, we can express our extended propagator as the 
upper left corner of an inverse matrix: 



i 



u>fj, — H 



(32) 



-i rs,r s 



By simple matrix partitioning we now find: 
-l 



-H, 



uji-n -n u „, .. 

=bb = bb 



(33) 
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This equation already reveals the algebraic structure we were aiming at. Of course, 
we can make it look more like Dyson's equation by introducing the matrix 



(34) 



which is easily recognized as the zeroth order of "H apart from infinitesimals 0(rf). 
When we introduce the "self energy" for the extended propagator by 

S(u) = H -e + H , 1 —7^K (35) 

=V =aa = =ab UJfi _ ft =ba 

=66 =66 

Eq. (p3[) becomes Dyson's equation for the extended two-particle propagator: 



wl-e- S(u) (36) 



Obviously, the zeroth order of the extended propagator takes on the following ap- 
pearance: 

£ (0) M]~ 1= <4-I (37) 

where the infinitesimal irj(2n r n s — 1) has been left out for simplicity. The anti- 
occupation number n r takes on the value 1 if the index r refers to an orbital unoc- 
cupied in the zeroth order ground state (virtual orbital) and if otherwise. 

Finally we can express the Dyson equation (|3^) in the appealing equivalent form: 

g(w) - g (0) (w) + g (0) (w)£(w)gM (38) 

We note that the self energy S_(cu) can be calculated from Eq. (|35|). Of course, 
the choice of the extended basis is not unique. Different ways of constructing finite 
subsets of this basis will hence lead to different approximations for the self energy. 
For example, powerful approximations may be achieved by constructing the basis 
set from correlated excited states, classified by the nature of the excitation with 
respect to the ground state. This is a common technique used in many modern 
many-body methods [ p5|]2(| . 

By construction, the extended Green's function 0_(uj) contains the full particle- 
particle propagator and, therefore, it can be used not only for scattering but also for 
calculating double ionization energies and the corresponding transition amplitudes. 
One has to be careful, however, when employing perturbation theory to extract 
approximations to the self energy S_(u>) for this purpose. The reason is that there 
are states in the extended basis contributing to Tt,, and degenerate in zeroth order 

° ^66 ° 

with states \Y rs ) where r and s are hole indices (i. e., assign occupied orbitals in the 
unperturbed ground state). It is, of course, possible to obtain valid approximations 
by applying the formalism of degenerate perturbation analysis. The deeper reason 
for the appearance of these degeneracies lies in the indefinite choice of the metric /t 
in Eq. (|20|). It is an open question, whether suitable extensions \rs m ( P^| ) can be 
found that avoid the degeneracies by using a definite metric. For a more detailed 
discussion of the properties of GioS), of perturbation theory, and of different choices 
of extensions Xrs see reference |24] . 
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D. Resolvents and Energy Dependence 



A striking new property of the extended propagator is revealed when we com- 
pare the various one- and two-particle propagators and try to interpret them as 
resolvents of some operator. We remark that, of course, all the matrices in the two- 
particle quantum numbers (r, s) are not matrices in a literal sense, since they are 
in general of infinite dimension and the indices may be continuous. These matrices 
can rather be understood as linear operators in the two particle Hilbert space. The 
matrix e is then a representation of the zeroth order Hamiltonian H D , restricted 

to the two particle space. In this sense, the zeroth order propagator Q^°'(ui) is 
a resolvent of this Hamiltonian, as Eq. J37| ) shows. Of course, the full extended 
propagator Q_{lo) can be seen as a resolvent of the energy-dependent operator 

g + 

In section HIE we will see that indeed this operator is an effective two particle 
Hamiltonian that can be used to describe elastic two-particle scattering off a cor- 
related 7V-fermion target. We want to remind the reader that neither the particle- 
particle propagator Tl(u>), nor one of its components n ± (w) is a resol vent in the 



complete two particle space, like the discussion at the beginning of [I B shows. 
Still, in a restricted space these propagators may be understood as resolvents. But 
the energy-dependent effective Hamiltonians that could be derived this way are 
pathological, since they contain linear contributions in lo that diverge at high en- 
ergies. The energy dependence of the two-particle self energy S_(uj), on the other 
hand, is given by Eq. (j35[), where the analytic structure is determined through a 
partial fraction series. For practical purposes this implies that the self energy has 
a finite limit for values of the energy variable oj far outside the spectrum of the 
Hamiltonian. 



E. One— Particle Potentials 



Another new feature of the extended two-particle propagators concerns the be- 
haviour with respect to one-particle potentials. For a system interacting by one- 
particle potentials the Hamiltonian H is a one-particle operator, i. e. it can be 
written in the form 

H = Y^ h n a \ a j ( 39 ) 

ij 

In this case, the particles of the system do not interact by real two-body forces 
but only with external or mean field potentials. Any operator of the form ([39]) can 
be diagonalized (i. e. hij oc Sij) by a simple unitary transformation Uj of the set 
of orbitals (basis of one-particle wavefunctions) . In general, this transformation 
includes the full one particle space. It can be seen directly from the definitions 
(pi|), (|l|), jl^), (|l4|), and (plf), that all of the propagators we have introduced so far, 
transform with one unitary transformation matrix on each one-particle index. For 
example, the particle-particle propagator transforms like: 

flij,kl = O tjj> N-i'j' ,k'V tk'kU'l (40) 

i' j' k'V 

Since a diagonalized one particle Hamiltonian is equivalent to the zeroth order 
Hamiltonian (i. e. it can be written in the form of Eq. @), we can easily study the 
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behaviour of the various propagators under the influence of a one particle poten- 
tial by performing the inverse transformation on the zeroth order propagators. It 
is particularly interesting to observe what happens under this procedure to equa- 
tions J3^ ) and ([36]) for our extended two-particle propagator, and Eq. (0) for the 
single-particle Green's function. From the argumentation given above, it follows 
immediately that in both cases the self energies S_(lu) and become energy in- 

dependent for one-particle potentials. It is easy to see that they are exact already 
in first order in the one-body interaction. For completeness we give £(w) and S_(uj) 
explicitly in these cases: 

SpqM = h pq (41) 
&rs,r' s' (^-0 — h rr / S ss ' ~{~ h ss ' S rr ' h rs i 8 ST ' h sr i S rs ' (^^) 

These matrices can be seen as representations of the one-particle Hamiltonian in 
the spaces spanned by one-particle wave functions aj |vac) and two-particle Slater 

determinants a\,,a\, |vac), respectively. 

Speaking in terms of resolvents, we see that both our extended two-particle propa- 
gator as well as the single-particle Green's function now become resolvents of energy 
independent effective Hamiltonians and these Hamiltonians contain no higher than 
first order contributions in the interaction. We want to stress that this is not the case 
for the particle-particle propagator TI{lo) from ( p"3| ) or any of its constituents (uj). 
We already pointed out that these propagators can be viewed as resolvents only in 
restricted index spaces. Therefore, the behaviour of their matrix inverses under 
transformations of the orbital basis is complicated and leads to energy-dependent 
effective Hamiltonians with contributions in all orders of perturbation theory. For 
this reason the calculation of scattering amplitudes or double ionization energies via 
the particle-particle propagator or its constituents in the presence of one-particle 
potentials is formally nearly as complicated as the treatment of a fully correlated 
system (as long as one does not transform to new orbitals, of course). When us- 
ing our new extended two-particle propagators, however, already the first order 
contribution to its self energy yields the influence of the one-particle potentials 
exactly. 



F. Time Space Equations 

For the sake of conceptual clearness we have developed the algebra of the extended 
two-particle propagators in energy space. In order to apply our developments to 



the scattering problem of section III , we will consider briefly the Fourier transforms 



into time space. We use the following convention for Fourier transformation: 

/OO 7 
^ e " iWt /H ( 43 ) 

In order to agree with standard notation we define: 

Qrs,r's' (^5 ^ ) ■ — Qrs.r' s' (t ^ ) 
Ilrs^'s' (^j t ) .= Yl rSjr ' s ' (t t ) 

When transforming the extended two particle propagator £/ rSir v(w), it is most 
convenient to start from Eq. ( |26| ) and make use of the convolution theorem. The 
result reads: 



10 



Qrs.r' s' {pi^ ) n^s^'g/ (tj t ) 

+&(t'-t) | i(vac| a s (t)4(t') |vac) • aj,(i')a r (t) |* )| 
-{r < — > s} - {r' < — > s'} + {r < — ► s , r' < — ► s'} (44) 



where the Fourier transform of the particle-particle propagator from Eq. ( p. 3D is 
given by: 

n ra , r v(t,t / ) = -6(t-t')(V \ a r (t)a s (t)al,(t')al,{t') \V ) 
i 

-~0(f - t)<* | al,(t')al(t')a r (t)a s (t) (45) 

The two terms on the right hand side of the above equation represent the parts 
Ilj. , , (t, t') and n~ , , (t, t'), respectively. In these expressions the operators ap- 
pear in the Heisenberg picture with a time dependence given by 

i(t) = e i6t Ae~ lilt . (46) 
Heavyside's theta function is defined as follows: 

ntj\ f 1 for t > 

^) : ={ fort<0 < 47 > 

As usual, the signs of the infinitesimals ±irj bring about the time ordering via theta 
functions. Although the time ordering does have relevance in scattering by giving 
the right boundary conditions, it does not affect the algebraic properties, since 
Dyson's equation transforms, regardless of the chosen time ordering of Q , into 

^, r v(M') = ft 5 -(M') 

/OO pOO 
dt x I dt 2 g(% q {tM)S Pq , P > q >{tiM)Gp- q -,r-s>{t'i,t l ) (48) 
- OO J — OO 

p* q' 

The zeroth order propagator in time space reads: 

■[n r n s 0{t - t') - (1 - n r n s )6(t' - t)] (49) 



III. OPTICAL MODEL DESCRIPTION OF TWO-PARTICLE 

SCATTERING 

In this chapter we will show how an extended two-particle propagator that satis- 
fies Dyson's equation can be used to derive effective two-particle equations for the 
elastic scattering of a two-fermion projectile by a correlated many-body target. 
The extended propagator introduced in the preceding chapter serves as one exam- 
ple, but other two particle Green's functions can be used as well, provided they fulfil 
certain requirements which will be discussed. We will show how the improper self 
energy of this propagator can be identified with the scattering T-matrix. Further- 
more, we will derive an effective Lippmann-Schwinger equation and an eigenvalue 
equation involving a two-particle scattering amplitude and an optical potential that 
is given by the two-particle self energy. 
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A. The Scattering Problem 



The type of scattering process that can be straightforwardly described by the 
considered type of two particle propagators is the elastic scattering of a two particle 
projectile off a correlated many-body target. For simplicity we assume that all 
particles are indistinguishable fermions, although extensions to composite particle 
systems, e. g. a positronium projectile, are easy to do and usually simplify the 
actual calculations. The target is taken to be in its nondegenerate iV-fermion 
ground state \^ ) with energy . In the simplest case, the projectile may consist 
of two (asymptotically) free fermions and is denoted by: 

a p a) q | vac) 

It is an eigenstate of the 'free particle' Hamiltonian which can be identified with the 
zeroth order Hamiltonian H from Eq. (||). All relevant information on the scatter- 
ing process is now contained in the two particle part n + (aj) of the particle-particle 
propagator from Eq. ([Ti]). However, we may consider a more general projectile 
consisting of an arbitrary two particle eigenstate \K) of some "free projectile" 
Hamiltonian H p . 

H p \K) = £ K \K) (50) 
\K) := A^ K |vac) (51) 

with A' K denoting the projectile creation operator: 

A< = J dpdq 4> K (p, q) aJ4 ( 52 ) 

Here the cumulative index K denotes a complete set of quantum numbers that 
characterize the two particle state \K) with energy Sr. Let p and q denote one- 
particle quantum numbers, e. g. momentum and spin p — (p,a). Since it appears 
more natural to use continuous sets of quantum numbers rather than discrete orbital 
indices when describing scattering processes, we changed the notation and replaced 
sums over one or two particle indices by integrals. This is merely a formal step, of 
course, since the integrals still contain discrete sums, e. g. over spin indices. The 
operators a p and a q create one particle states with the given quantum numbers. 
The quantity (f>K (p, q) plays the role of a two particle wavefunction. 

The simple case of a projectile consisting of two free particles now emerges as 
a special case in the more general treatment. The free projectile Hamiltonian H p 
reduces to the kinetic energy operator and we can take a p and a q to create plane 
wave states. In that case we can choose the wavefunction to be of the form 



<t>K{p,q) = 5 klP 8k 2 q (53) 

where K = (ki,k,2). In the following we will develop the theory for a general 
projectile state, but come back to the simpler case at times. 

The principal quantity for describing the results of scattering experiments is the 
so called scattering matrix or S-matrix: 

S(K ^ K>) = (54) 

We want to describe the elastic scattering of two particle states \K). The S-matrix 
then gives the probability amplitude for the scattering of the two particle projectile 
\K') off the target state \^f ), into the projectile state \K), the target remaining 
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in the state |\& ). Formula (54) gives the S-matrix as the overlap of the so called 
scattering states: 

|*|> = lim AUt)\*o)e- i£Kt (55) 

t — >^00 

Here the projectile creation operator A^ K appears in the Heisenberg picture defined 
in Eq. (ffq). The scattering states can be understood as states in the Heisenberg 
picture that behave asymptotically in the distant past/future like a free projectile 
\K) and target 

Strictly speaking, the limit in fl55| ) does not converge in the usual sense and is 
meant as an Abelian limit (27) : 

/•Too 

|*±) = nm/ dT e ±eT A ] k {t) \V )e~ l£KT (56) 

In the following we will use the Abelian limit where appropriate, without referring 
explicitly to it. 

For the S-matrix we find: 

S(K <- K') = lim e i£Kt {^ \ A K (t)A^ K ,(t') \^ }e- i£ ^ t ' (57) 



Employing the definition ( p2[ ) of the projectile creation operator A^ K , we arrive at 
the following expression: 

S(K<-K')= lim e l(£Kt - £ K' t,] J dpdqdp 1 dq' ■ 

■<j>* K (p,q)$> \ a q {t)a P {t)al,{t')a\,{t') \V ) <f> K ,(p',q') (58) 
= lim e ^Kt-e K ,t') I dpdqdp > dq < . 

t' — ► — oa ** 

■<t>K{p,q)^ plq ,(t,t')<p K ,{p',q') (59) 

The last equality follows the definition (^). This shows that the two particle part 
n + (£,i') of the particle-particle propagator contains all relevant information on the 
scattering process. Obviously, U_ + (t,t r ) can be replaced by G_{t,t') as well as by 
n(t, t'), since the additional terms do not contribute due to the fixed time ordering 
in equation (p^, like it can be seen from equations (H) and (45). 



B. General requirements for a scattering propagator 

The main objective of this chapter is to show how the extended two particle 
propagator Q_ of Eq. ( p4| ) or (|44|) can be used to reduce the many-body scattering 
problem to effective two-particle scattering equations. The derivation we will give 
in the following, however, only makes use of some fundamental properties of the 
extended propagator. Loosely speaking these are: the propagator used contains n + 
and satisfies Dyson's equation. Since it is possible to find different extensions of the 
particle-particle propagator that meet these requirements, we will introduce the new 
symbol Qk,k' {t, t') for a general propagator that fulfils certain requirements stated 
below. We will refer to the general two particle Green's function GK,K'(t,t') as the 
scattering propagator. For convenience this propagator is labelled with two particle 
quantum numbers like K and K' ', each of these replacing two orbital indices. The 
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transformation between these two methods of labelling two particle propagators is 
given by the two particle wavefunctions <j>x (p, q) like in Eq. ( fi2| ) . In the following 
we will state the required properties explicitly: 

The first requirement allows one to rewrite Eq. ( |57| ) for the S-matrix in terms of 
the scattering propagator: 

(i) For t > t' the scattering propagator reduces to the transformed two particle 
propagator U_ + (t, t'): 

iG K ,K'{t,t') = <*„| A K (t)A^ K ,(t') for t > t' (60) 

The second property requires the main new feature of our extended two particle 
propagator: Dyson's equation. It introduces the self energy S of the scattering 
propagator Q, which will serve as the two particle optical potential: 

(ii) The scattering propagator satisfies a Dyson equation: 

gK,K>(t,t>) = g% ) K ,(t,t>) 

+ J dKidK 2 dt l dt2g^ Kl {t,tl)S Kl ,K 2 (tl,t2)GK 2 ,K'(t2,t') (61) 

Here, a free projectile propagator Q^l appears that reads 

iSjK,K'M = *(* - th > 6K > K ' (62) 

for K or K 1 belonging to the subset of two particle quantum numbers used 
to describe the projectile states. 

The integrals over Ki in (|6l]) are meant to run over the complete two particle space. 
The 6 function denotes a proper normalization condition. The actual form of this 
condition depends on the set of quantum numbers used. In the example of two free 
projectile particles this normalization condition reads: 

$k,k> — ^fcife'j Sk 2 k' 2 - <5feife^ 8k 2 k[ (63) 

The second term on the right hand side accounts for the fermionic nature of the 
projectile's components. In the case of a bound state projectile, on the other hand, 
the index K may combine an internal quantum number n of the projectile with a 
center of mass wavevector k and consequently the S function is given by: 

5k, k< = $ nn < (64) 

The third requirement is a minor technical point. It is a common feature of 
many-body Green's functions depending on two time variables that these functions 
actually depend only on the time difference, when the total energy of the given 
system is conserved. Thus we require: 

(iii) As a consequence of the conservation of energy we want the function Q and 
the other functions that appear in the Dyson equation to depend only on the 
time difference t — t' . 

It is easy to see that these conditions are fulfilled by the extended propagator 
Gk 1 k 2 ,k' 1 k' 2 {t,t') 01 Eq. (fti|), in the case of a projectile consisting of two (asymp- 
totically) free particles where K = (fci,^) holds. The zeroth order propagator 
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^jfcil 2 k' k' *') f rom Eq. ( ^9| ) now can be identified with the free projectile propa- 
gator G^ K ,(t,f). 

In the more general case of a projectile in an arbitrary, possibly bound, two- 
particle state, the transformed function 

G K ,K>(t,t') := J dpdqdp' dq' <t > * K {p,q)g pq . p , q ,(t,t')<f >K ,(p',q') (65) 

fulfils the requirements. We only have to be a little more careful, if the projectile 
state \K) is not a Slater determinant and thus cannot be expressed as the eigenstate 
of a one-particle operator. This is generally true for bound two-particle projectile 
states. Since the zeroth order Hamiltonian H Q of Eq. (||) is usually required to be a 
one-particle operator, the free projectile propagator of Eq. (|6^ ) will generally 
differ from the zeroth order of Eq. ( |4S| ) both in the phases and the set of 
quantum numbers that diagonalizes the propagator. This does not cause a serious 



problem, however, since the Dyson equation (48) for Q can be rewritten with 
replaced by the propagator . This implies a trivial renormalization of the u> 
independent part of the self energy S_(co), as the energy-space equations (am and 



( |38[ ) show. 

If we still want to keep the direct connection to perturbation theory, such that 
is the zeroth order of Q, we have to find some way to include H p , which is 
primarily an operator in the two particle Hilbert space, into the zeroth order (Fock- 
space) Hamiltonian H Q . In this case we have to take care that the free projectile 
states do not interfere with the zeroth order ground state of the target, which also 
has to be an eigenstate of H . This could be done by constructing the projectile 
from unoc cup ied (virtual) orbitals only. A procedure similar to this is described in 
reference fl9|| . 

We now carry on without referring to the extended propagator introduced in the 
last chapter but assume that we have a Green's function, the scattering propagator 
G, which possesses the properties stated above. 



C. Relation to the Scattering T— Matrix 

Now we are in the position to link up with Eq. j5^ ) and express the S-matrix in 
terms of the scattering propagator using property (i): 

S{K^K')= lim jfe, K -(M')e !(£K '^ K ' t ' ) (66) 

In order to see the connection to the popular T-matrix description of scattering 
we rewrite Dyson's equation (|6l| ) introducing the improper self energy or Green's 
function T-matrix TK,K'{t, t')- 

g = g ip) + g {p) fg ip) 

= = .7 b) r = (67) 
T = £ + £g w T 

Here, we are using a shorthand notation where the entities like Q_ are matrices with 
indices (K,t) and (K',t'). The matrix product then implies integration over all 
inner indices, namely the two particle quantum numbers and the time variable. For 
example: 
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(K,t)(K>,t>) 



J ^i^rfti^^^^ti)^,^^,^)^^^,*') (68) 

It is easily found that inserting ( |67| ) into the expression (|6^) for the S-matrix and 
using property (i) yields: 

S(K K') = &kk' - 2mf K , K ,{£ K ) 5(£ K - £k>) (69) 

Here, the quantity Tk,K'{£k) denotes the Fourier transform of the improper self 
energy TK,K'(t,t'), evaluated at energy £k- Transforming T_ on a single energy 
variable becomes possible, because the improper self energy is a function of the 
time difference only. This is a consequence of property (iii) . 

Eq. ( |69|) now allows us to identify the improper self energy 7k } k'(£k) with the 
on-shell scattering T-matrix t(K <— K'), which is defined by 

S(K <- K') = 8 K k> - 2m 8{£ K - £ K ,) t{K <- K'). (70) 

This shows a connection between the many-body Green's function Q and effective 
scattering quantities represented by the common scattering S- and T-matrices, 
which is analogous to the problem of scattering one fermion off a correlated many 
fermion target p8 29 1. In that case the single-particle Green's function G(t,t') 



takes on the role of Q(t, t'). 



D. Effective Lippmann— Schwinger Equation 

We will now proceed to further carry this analogy and derive effective two-particle 
scattering equations. First we define scattering amplitudes that take the role of 
effective two-particle wavefunctions: 

f K '+(K,t):=(* \A K (t)\*+,) (71) 

These amplitudes contain all the relevant information concerning the scattering 
process, since we can express the S-matrix through: 

S(K <- K') = lim e l£Kt f K ' + (K,t) (72) 

t — ^OO 

The relation to the Green's function is given by: 

f K '+(K,t)= lim iG^K'Me-^'*' (73) 

t' — > — oo 

The time dependence of the scattering amplitudes can be taken from the defini- 
tion, using the fact that |^ ) and are eigenf unctions of the Hamiltonian: 

pK'+flf +\ — AT; I A.,l+\ UT/ + 



f" + {K,t) = (p \A K (f)\*+,) 



= {y o \e im A K e- lHt \*j 



- -iS K ,tfK'+ 



e 



K*f" +(K) (74) 



The last equality defines the time independent scattering amplitude. We note that 

f K '+(K)=e i£ «' t f K ' + (K,t) (75) 
= (* \A K \*+,) (76) 
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We can now use Dyson's equation for Q to derive an effective Lippmann- 
Schwinger equation for the time independent scattering amplitude / + (K). In- 
serting Dyson's equation ( |6l| ) into ([73]) and using ([75|) leads after some simple 
manipulations to the equations: 

f K ' + {K)= lim ie if *''S^(t,i')e- if *''' 

V — > — oo 

= $kk' + lim i / dKi dK 2 dt\ dt2 x 

i'-»-oo J 

e i£ x ,t g(p) Ki (tj tl) 5 Xi)jKa ( tli t2 ) g K2iK ,(t 2 ,t') e-**'* (77) 

Carrying out the time integrations with the help of property (iii) leads to the desired 
working equation for the time independent scattering amplitude: 

f K ' + (K) = 6 K k> + J dK 1 dK 2 g^ ) Ki {£ KI )S K ^ K2 {£ KI )f K ' + {K 2 ) (78) 

This is an integral equation which can be solved for the scattering amplitude once 
the self energy S_(lo) is known. It is analogous to the Lippmann-Schwinger equation 
of basic scattering theory. In the above equation (78) the self energy replaces the 
potential that appears in the Lippmann-Schwinger equation. Therefore we have 
actually derived an exact optical model description for our many-body scattering 
problem: The initial N + 2 particle problem has been reduced to the problem of 
scattering two particles by an optical potential. The optical potential, which is 
given by the self energy S_(lo), is in general energy dependent and takes account of 
all the correlation within the target as well as the interaction between the target and 



the projectile particles. Like we have seen in section HE, this self energy becomes 
energy independent for systems with particles interacting by one particle potentials 
only. 

All this is analogous to the elastic scattering of one fermion by a many-body 
target. In this case the self energy S(w) belonging to the single-particle Green's 



function G(u>) takes the role of the optical potential [E8 29 q,M], which becomes 



energy independent for systems interacting by a one particle potential. 



E. A Two— Particle Eigenvalue Equation 

In basic scattering theory the Lippmann-Schwinger integral equation presents an 
alternative formulation to Schrodinger's eigenvalue equation. Thus we can expect 
that in an optical model description it is, in general, possible to formulate an ef- 
fective eigenvalue equation for the scattering amplitudes. Multiplying the integral 
equation ([78]) with the energy difference (£k> — £k) and using the Fourier transform 
of ( |62|) leads us to the desired equation: 

(£k -£ K >)f K ' + {K) + J dK 2 S K , K2 {£ K ,)f K ' + {K 2 ) = (79) 

This is an eigenvalue equation in the two particle Hilbert space. It is formulated us- 
ing a basis of two particle wavefunctions that diagonalize the projectile Hamiltonian 
Hp defined in III A| . Reintroducing H p and switching from the quantum number 



representation to abstract two-particle-space vectors we can rewrite ( |78| ) into: 

[H p + S(£ K ,)} \f K ' + ) = £ K , \f K '+) (80) 



17 



We have defined the effective scattering kets \f K +) with the help of the free pro- 
jectile states \K) of Eq. (|5l| ) by 

(K\f K '+) = f K ' + (K) (81) 

The operator S(Sk') is defined by it's matrix elements: 

(Ki\S(w)\K 2 )=S KltK ,(w) (82) 

Of course, the Green's function Gk 1 ,k 2 {< jJ ) is now rediscovered as a representation of 
a resolvent Q{oS) that satisfies the following operator identity in two particle Hilbert 
space: 

{h p + S(u) - u\ g(u) = id (83) 



This links up with section 
ticle propagator as a reso 



II D . where we already interpreted the extended two par- 
vent of an uj dependent operator. Now we have finally 
proven that the scattering propagator (and therefore also the extended two particle 
propagator from the last chapter) actually is the resolvent of an effective Hamilto- 



nian that allows to cal culat e S-matrix elements through (pOJ) and (72)! Naturally 



the remarks of section II D on the energy dependence of the effective Hamiltonian 
apply here in the same way. 

We want to remind of the fact that the effective Hamiltonian and therefore also 
the optical potential derived directly by Feshbach projection contains a divergence 
that is avoided by our two particle self ener gy. A lso the behaviour with respect 



to one particle potentials discussed in section HE provides a contrast between the 



two approaches. This leads us to the following conclusion: Given some arbit rary 



two particle propagator, it is the properties (i), (ii), and (iii) of section HI B that 
lead to the effective scattering equations of this last chapter. The usefulness of 
these equations for approximative calculations of scattering parameters, however, 
is determined by the properties of its self energy. As shown in chapter |TJ, the 
self energy related to our extended propagator possesses a good-natured energy 
dependence and behaves well with respect to one particle potentials. For these 
reasons the extended two particle propagator provides a useful tool in scattering. 
Furthermore it presents an alternative approach that removes drawbacks of the 
Feshbach projection effective potentials discussed in reference p9|. 



IV. CONCLUSIONS 



In this paper we have presented an extension of the fermionic particle-particle 
propagator that has algebraic properties similar to the single-particle Green's func- 
tion. We explicitly constructed one example of an extended propagator and in- 
vestigated its major algebraic properties. These include an analogon of Dyson's 
equation with a two-particle self energy We showed that this self energy serves as 
an exact optical potential for the scattering of a two-particle projectile by a corre- 
lated many-body target. The projectile may either consist of two asymptotically 
free fermions or a bound two-particle state. In the latter case the two-particle self 
energy has to be trivially renormalized. 

In contrast to effective two particle potentials derived directly from the multi- 
ple particle Schrodinger equation, the two particle self energy has a good-natured 
energy dependence. Additionally, correlations due to a one-particle potential are 
accounted for exactly already by a first order, energy-independent self energy. 
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The ideas behind the construction of extended two-particle propagators are quite 
general and lead to a whole new class of extended two-particle propagators with 
similar algebraic properties. Employing approximations that aim at ab initio cal- 
culations will lead to different working equations for each of these functions. A new 
feature of the extended propagators is the appearance of degeneracies that may 
require the use of degenerate techniques necessary when employing perturbation 
theory based approximations. The reason for this lies in the indcfmitcncss of the 
metric used in the construction of the extended propagators. It is an open ques- 
tion, whether or not an extended propagator can be found that combines a definite 
metric with the advantageous properties we discussed in the present paper. 

The present approach can be generalized to construct optical potentials for the 
scattering of an M-particle projectile by a many-body target. The starting point 
is the introduction of an extended M-particle propagator which fulfills a Dyson 
equation. 

NOTES 

1. We use the following definition for the scalar product, which is canonical: 

(('«)'('« ))-<*w + w« 

The action of a diagonal operator is defined by: 

[A, \ / >= ( A 1 \<j> l )\ 



2. Instead of the matrix partitioning carried out in |2l| , one can equally well 
proceed by introducing projection operators onto the subspace spanned by 
the set of vectors { |V" 9 G )} and its orthogonal complement pQ[ . 

3. For these states the rules of footnote [l] apply likewise. Additionally we define 
for product states: 

(\a')®{b'\, \a) ® {b\ > = {a'\a} ■ (b\b') 
and (A®B)(\a)®{b\)=A\a)®(b\B 

4. For earlier work on effective two-particle scattering equations in a different 
context see p9| . 
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